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Abstract 


We study the P — V criticality and phase transition in the extended phase space 
of charged Gauss-Bonnet black holes in anti-de Sitter space, where the cosmological 
constant appears as a dynamical pressure of the system and its conjugate quantity 
is the thermodynamic volume of the black holes. The black holes can have a Ricci 
flat (k = 0), spherical (k = 1), or hyperbolic (k = —1) horizon. We find that for the 
Ricci flat and hyperbolic Gauss-Bonnet black holes, no P — V criticality and phase 
transition appear, while for the black holes with a spherical horizon, even when the 
charge of the black hole is absent, the P—V criticality and the small black hole/large 
black hole phase transition will appear, but it happens only in d = 5 dimensions; 
when the charge does not vanish, the P—V criticality and the small black hole/large 
phase transition always appear in d = 5 dimensions; in the case of d > 6, to have the 
P — V criticality and the small black hole/large black hole phase transition, there 
exists an upper bound for the parameter b = a|Q|-2/ (4-3) where & is the Gauss- 
Bonnet coefficient and Q is the charge of the black hole. We calculate the critical 
exponents at the critical point and find that for all cases, they are the same as those 
in the van der Waals liquid-gas system. 
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1 Introduction 


Over the past years, a lot of attention has been attracted to study thermodynamic prop- 
erties of black holes in anti-de Sitter (AdS) space mainly due to the AdS/CFT corre- 
spondence |1, 2, 3], which says that thermodynamics of black holes in AdS space can be 
identified with that of dual strongly coupled conformal field theory (CFT) in the bound- 
ary of the AdS space. The thermodynamic properties of black holes in AdS space are 
quite different from those of black holes in asymptotically flat spacetime or in de Sitter 
space. In AdS space, large black holes are thermodynamically stable, while small ones are 
thermodynamically unstable. Under a certain temperature, there does not exist any black 
hole solution in AdS space, and a so-called Hawking-Page phase transition can happen 
between stable large black holes and thermal gas in AdS space [4]. This phase transition 
can be explained as the confinement/deconfinement phase transition of gauge field in the 
AdS/CFT correspondence [5] 

Thermodynamics of charged black holes in AdS space has also been intensively studied 
in the literature. In particular, in a canonical ensemble with a fixed charge of black holes, 
an interesting phase transition between large black hole and small black hole has been 
found in [6]. A critical point exists in the phase diagram with charge Q and temperature 
T of the black holes. The phase transition and critical behavior of the charged black holes 
are quite similar to those in the van der Waals liquid-gas system. Note that although the 
Q- diagram of the charged black holes is similar to the P-V diagram of the van der Waals 
system [7, 8], where ® is the chemical potential conjugate to the charge Q of the black 
hole, while P and V are the pressure and volume of the van der Waals system, the analogy 
is problematic since the charge Q is an extensive quantity and ® is an intensive one in 
black hole thermodynamics, while P is an intensive quantity and V is an extensive one in 
the van der Waals system. 


Recently the analogy of charged black hole in AdS space as a van der Waals system has 
been further enhanced by studying critical behavior of the charged AdS black hole in the 
extended phase space [8], where the cosmological constant appears as the thermodynamic 
pressure and its conjugate quantity as a thermodynamic volume of the balck hole. It 
has been showed that both systems share the same critical exponents and have extremely 
similar phase diagrams. This analogy has been generalized to the higher dimensional 
charged black holes, rotating black holes and Born-Infeld black holes in AdS space [9]. See 
also relevant discussions in various cases in the references [10]. For a most recent work 
with further interesting observations on the higher dimensional rotating AdS black hole 
with a single rotation parameter, see [11]. 

In the usual discussions of thermodynamic properties of black holes in (A)dS spaces, 
the cosmological constant is treated as a fixed parameter. Indeed different cosmological 
constant implies a different gravity theory under consideration. In the first law of black 
hole thermodynamics, the mass M, angular momentum J and charge Q are conserved 
charges, while the cosmological constant is a parameter for a given theory. However, there 
are some physical reasons to view the cosmological constant as a variable [8]. For example, 
one may suppose that there exist some more fundamental theories where some physical 
constants such as Yukawa coupling, gauge coupling constants, Newton’s constant, and/or 
cosmological constant, may not be fixed values but dynamical ones arising from the vacuum 
expectation values. In that case, it is natural to add these variations of “constants” into the 
first law of black hole thermodynamics [12, 13]. For example, in the gauged supergravity, 
the cosmological constant appears as a coupling constant. In addition, if the cosmological 
constant term is absent in the first law of black hole thermodynamics, the scaling argument 
cannot lead to a consistent Smarr relation for balck hole thermodynamics. Similar situation 
appears for the Born-Infeld black holes. To get a consistent Smarr relation by scaling 
argument, one has to introduce the Born-Infeld parameter term into the first law of Born- 
Infeld black holes [14, 15]. Once one regards the cosmological constant as thermodynamic 
pressure in the first law, the black hole mass M should be explained as enthalpy rather 
than internal energy of the system [16]. In the geometric units Gy = h=c=k=1, one 
can identify the cosmological constant with the pressure as 


i, _ -p= 


a lén2z (1) 


in d dimensional spacetime. Then the thermodynamic quantity conjugate to the pressure 
is called “thermodynamic volume” of black holes [17, 18, 19, 20, 21]. 

With the new identification, one can study more thermodynamic quantities of black 
holes such as adiabatic compressibility, specific heat at constant pressure, or even the 
“speed of sound” associated with the black holes [22, 17, 23]. Furthermore, as pointed 
out by Dolan [17], this also leads to an interesting possibility to reconsider the critical 
behavior of AdS black holes in an extended phase space, including pressure and volume 
as thermodynamic variables. In this way, one can more naturally consider the analogy 
between the AdS black holes and the van der Waals liquid-gas system. Indeed, the authors 
of Ref. [8] initiated the investigation of the P — V critical behavior of a charged AdS black 
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hole in the extended phase space and found that the P — V diagram of the black hole is 
exactly the same as the one for the van der Waals liquid-gas system. The phase transition 
of the charged AdS black hole has the same critical exponents as the van der Waals system. 
As a result the analogy between the charged AdS black hole and the van der Waals system 
becomes complete. 

In this paper we are going to study the P — V criticality in the extended phase space 
of charged Gauss-Bonnet black holes in AdS spaces. The motivations are as follows. 1) 
Is there any effect of some higher derivative terms of curvature on the P — V criticality 
in the picture of [8, 9]? 2) Note that in the discussions of [8, 9, 10], the electric charge Q 
of the black holes plays an important role in the analogy, we want to know whether the 
charge is essential in this picture. 3) In the above discussions, the topology of the black 
hole horizon is a sphere. Note that in AdS space the black hole horizon could also be 
Ricci flat or hyperbolic [24]. Therefore it is interesting to see whether the P — V criticality 
appears or not for those topological black holes. Indeed we find some new features when 
the Gauss-Bonnet term is present. 

This paper is organized as follows. In the next section we give some thermodynamic 
quantities of the charged Gauss-Bonnet black holes in AdS space. In section 3, we will 
investigate the P — V criticality of the Gauss-Bonnet black holes without charge, and 
compare with the van der Waals system. In section 4, we will discuss the general case with 
charge. Section 5 is devoted to the conclusions and discussions. 


2 Thermodynamics of Gauss-Bonnet black holes in 
AdS space 


Consider the d-dimensional Einstein-Maxwell theory with a Gauss-Bonnet term and a 
cosmological constant A: 


5 = = d*z./—g|R — 2A + acp(Ru R — 4R R” + R?) —4rF,,F""), (2) 
where agp is the Gauss-Bonnet coefficient with dimension [length]? and the cosmological 
constant A = ees Fv is the Maxwell field strength defined as F),, = 0, A, — O, Ap 
with vector potential A,,. In the low energy effective action of heterotic string theory, agg 
is proportional to the inverse string tension with positive coefficient [25]. Thus in this paper 
we will consider the case with a positive Gauss-Bonnet coefficient, namely, agg > 0. In 
addition, let us mention here that the Gauss-Bonnet term, Rap = Ront —4Rı R” + 
R?, is a topological term in d = 4 dimensions and has no dynamics in this case. Therefore 
we will consider d > 5 in what follows. 

The action admits a static black hole solution with metric 


ds? = —f(r)dt? + f (r)dr? + r°hijda dat, (3) 


where h;;dx'dz/ represents the line element of a (d — 2)-dimensional maximal symmetric 
Einstein space with constant curvature (d — 2)(d — 3)k and volume Nx. Without loss of 
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the generality, one may take k = 1, 0 and —1, corresponding to the spherical, Ricci flat 
and hyperbolic topology of the black hole horizon, respectively. The metric function f is 
given by [25, 26, 27, 28] 
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where & = (d — 3)(d — 4)agp, M is the black hole mass, Q is related to the charge of the 
black hole and P = -Ż. Note that in order to have a well-defined vacuum solution with 
M = Q = 0, the effective Gauss-Bonnet coefficient a and pressure P have to satisfy the 


following constraint 
64raP 
0 < —— <l. 5 
~ (d—1)(d—2) 7 (5) 
The horizon radius r, of the black hole is determined by the largest real root of the equation 
f (rn) = 0 and the mass M can be expressed in terms of the horizon rp, 
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The Hawking temperature of the black hole can be easily obtained by requiring the absence 
of conical singularity at the horizon in the Euclidean sector of the black hole solution, which 
is given by 

16nPré/(d — 2) + (d —3)kr? + (d —5) ka — 22, 
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Since we are going to discuss the thermodynamics of the black hole in the extended phase 
space by introducing the pressure P = -E, the black hole mass M should be considered as 
the enthalpy H = M rather than the internal energy of the gravitational system [16]. Other 
thermodynamic quantities can be obtained through thermodynamic identities. For exam- 
ples, the entropy S, thermodynamic volume V and electric potential (chemical potential) 
® are given by 
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We can see that the thermodynamic volume is a monotonic function of the horizon radius 
ra in our case. Taking advantage of this, one can use the horizon radius to stand for the 


conjugate volume in the equation of state. In the later discussions, we will see that just 
the horizon radius corresponds to the specific volume in the Van der Waals equation rather 
than the thermodynamic volume. 

It is easy to check that those thermodynamic quantities satisfy the following differential 
form 


dH = TdS + dQ + VdP + Ada, (11) 
“ee OH (d — 2)k?d (d — 2)kX,T 
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is the conjugate quantity to the Gauss-Bonnet coefficient &. Note that here we have treated 
the Gauss-Bonnet coefficient as a variable. By the scaling argument, we can obtain the 
generalized Smarr relation for the black holes 


(d—3)H = (d—2)TS —2PV + 2Aa + (d—3)Q®. (13) 


For black holes in Lovelock gravity, a generic relation is given in [31]. Note that for the 
Born-Infeld black hole case, the Born-Infeld parameter b has units of the electric intensity 
and the corresponding conjugate quantity $ has units of electric polarization, the authors 
of Ref. [9] therefore interpret the conjugate quantity B as Born-Infeld vacuum polarization. 
In our case, note that the Newton’s constant Gy has dimension [length]*~?, A in fact has 
dimension [length]~°. In addition, the black hole has volume V ~ Erit, area A ~ Erri ?, 
the black hole horizon has scalar curvature Rp ~ k/ r, and the Gauss-Bonnet term on the 
horizon Rep ~ k?°/rf. The first term in (12) then has the form V Rap, while the second 
term TAR. Clearly when k = 0, both terms vanish. In particular, it is quite interesting 
to note that the first term in (12) is just the second term in the black hole mass (6), while 
the second term in (12) is just the second term of the black hole entropy (8) multiplied by 
Hawking temperature T. This results is consistent with the one given in [31]. There the 
authors have shown that the conjugate quantity YU“) to the Lovelock coefficient by generally 
consists of three terms, which are related to mass, entropy and the anti-symmetric Killing- 
Lovelock potential of the black hole solution, respectively. Anyway, the conjugate quantity 
A is certainly related to the effective stress-energy tensor of the Gauss-Bonnet term Rep. 
It is required to further understand the physical meanings of the conjugate quantity A in 
the black hole thermodynamics in the extended phase space. 

Furthermore, the Gibbs free energy and Helmholtz free energy can be obtained by 
Legendre transformations as 


G=G(T,P,Q)=H-TS, F=F(T,V,Q)=G-— PV. (14) 


It is worth pointing out that the Helmholtz free energy F is calculated by subtracting 
off the contribution of the background of the AdS vacuum solution. Therefore, only the 
case with a negative F is regarded as the black hole solution being thermodynamically 
favored over the pure AdS vacuum solution. On the other hand, one can see that the 
black hole entropy (8) may be negative in the case with hyperbolic horizon (k = —1) [29]. 


A negative entropy does not make any sense in statistical physics. As a result, in the 
following discussions, we will impose the following constraints 
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In addition, let us note that in the case k = —1, one can see from the metric function f (4) 


that there exists a minimal horizon radius r? > 2. However, the non-negative definiteness 
of the black hole entropy (8) gives a more strong constraint on the horizon radius in this 
case: r? > (2+ 4/(d—4))a. 


3 P-V criticality in the case with Q = 0 


From Eq. (6) we can see that the charge Q controls the lowest power term of rp. So one 
might expect that the critical behavior would be different between the case Q = 0 and the 
case Q # 0. We will first study the case Q = 0 in this section. In this case, the equation 
(6) is reduced to 
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From the Hawking temperature (7) we have the equation of state of the black holes 
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Since we want to compare it with the van der Waals equation, we make a series expansion 
for the van der Waals equation with the inverse of specific volume v, 
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Thus one can identify the specific volume v with the horizon radius of the black holes 
as [9, 8] 
Ar), 
ale ea 
We see that the specific volume v is proportional to the horizon radius rp, therefore we 
will just use the horizon radius in the equation of state for the black holes hereafter in this 
paper. Of course, one can easily recover all the results in terms of v by the relation (19). 
We know that the critical point is determined as the inflection point in the P — V 
diagram, i.e., 


(19) 
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= 0, (20) 


where we have used the subscript ‘c’ to stand for the quantities at the critical point. From 
(17) we can obtain the critical temperature 
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and the equation for the critical horizon radius (specific volume) 
(d —3)r,,, — 12kr?.a@ + (12d — 60)a? = 0. (22) 


It can be seen that the physical solutions of equation (22) depend on the dimension and 
topology of the black hole horizon. We will give a detailed discussion for each situation in 
the next subsection. 


3.1 The phase transition for different horizon topologies and di- 
mensions 


The topology of the black hole horizon is characterized by the value of k, where k = 0, 
k = —1, and k = +1 correspond to the Ricci flat, hyperbolic and spherical horizon, 
respectively. For each case of k, we will investigate the phase structure and criticality in 
the extended phase space. 


1) Ricci flat case with k = 0 


For the Ricci flat case, i.e. k = 0, the equation (22) has no positive solution for rre, which 
means that there does not exist any critical point. This can be seen more clearly from the 
equation of state (17). When k = 0, it reduces to 

_d-2 


P= T. 2 
Ar), (23) 


The pressure as a function of r, is monotonic for a fixed temperature. Thus no phase 
transition can happen. 


2) hyperbolic case with k = —1 


In this case, the equation (22) can be rewritten as 
(d — 3)rp. + 12r} © + (12d — 60)a? = 0. (24) 


Only when a < 0 may this equation admit positive roots. Thus we conclude that there 
does not exist any phase transition and critical point in the hyperbolic case. We can also 
see the fact from the equation of state of the black holes 

d—2 2a d—2)(d—3 d — 2)(d—5)a 
_d=24 Šp, (@=2)(d=3) _ (d= 2)(d= 5) mm 
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Note that the fact r? > 2a, the first term in the right hand side of the above equation is 
always positive. When d = 5, the pressure is a monotonic function of rp, therefore there are 
no critical point and phase transition in this case. When d > 6, the last term in the right 
hand side of the above equation dominates for the small radius case. In this case, there is 
only one stable phase. Again, there do not exist any critical point and phase transition. 


3) spherical case with k = 1 
In this case, the equation (22) reduces to 
(d — 3)r;,, — 12r}. + (12d — 60)a? = 0. (26) 


This equation has a positive real solution for rac only in five dimensional and six dimensional 
cases. When d = 5, we have 


1 1 
Tho = V60, P.=—-, T,= = 27 
48TA TV 24a ( ) 
from which we can easily find an interesting relation among the critical pressure P., tem- 
perature T, and horizon radius rpe: 
Pj. 1 


a (28) 


This relation is universal in the sense that it is independent of a. This result is very similar 
to the one in the van der Waals system, where the critical point also has an analogue 
universal relation with the only difference that the right hand side of equation (28) is 3/8. 
From (28), we can also see that it is more natural to view the horizon radius as the specific 
volume, instead of the thermodynamic volume. The corresponding P — rp diagram drawn 
in the left plot of figure (1) is exactly the same as the P — V diagram of the van der Waals 
liquid-gas system. The critical isotherm is denoted by the red line in figure (1). For a fixed 
temperature lower than the critical one, we have two branches whose pressure decreases as 
the increase of horizon radius, one is in the small radius region (corresponding to fluid phase 
in the van der Waals system) and the other is in the large radius region (corresponding to 
the gas phase). Such two branches have a positive compression coefficient, thus representing 
stable phases. Between them there is an unstable phase with a negative compression 
coefficient. For appropriate values of pressure, the isothermal line allows two physical 
horizon radii. Therefore, the so-called small black hole/large black hole phase transition 
occurs, which is reminiscent of the liquid/gas phase transition of the van der Waals system. 
Such phase transition is first order for T < T,, while it becomes second order at T, just 
as the same as the case in the van der Waals system. Above the critical temperature, the 
black holes are always in the gas phase and no phase transition happens. 

In the d = 6 case, although we can get a positive solution of the horizon radius from 
(26), which has a corresponding positive critical temperature T,, there does not exist any 
phase transition below this critical temperature. As one can see from (17) for k = 1 and 
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d = 6, the last term now dominates the equation of state at small horizon radius, where 
the pressure increases with the radius. Thus, comparing to the d = 5 case, the stable phase 
with small horizon radius no longer exists. This can be seen clearly from the right plot of 
figure (1), which exhibits the P — r, diagram in the d = 6 case. We see from the plot that 
for any fixed temperature there is only one branch with a postive compression coefficient. 
In other words, there do not exist two stable phases between them a phase transition could 
occur in this case. 
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Figure 1: The P — V(r) diagram for d = 5 (left plot) and d = 6 (right plot) near the 
critical temperature, respectively. Without loss of generality, we consider the case a = 1. 
In both plots, the temperature of isotherms decreases from top to bottom and the pressure 
satisfies the relation (5). 


We plot the Gibbs free energy defined in (14) in figure (2) as a function of temperature 
for various pressures. We can see from the left plot that in the case of d = 5, the Gibbs 
free energy with respect to temperature develops a “swallow tail” for P < P., which is 
a typical feature in a first order phase transition, and at P = P., the “swallow tail” 
disappears, corresponding to the critical point. On the other hand, in the case of d = 6, 
there does not exist any “swallow tail” and thus no phase transition happens. 


3.2 Critical exponents at the critical point 


From the above discussions, even when the charge is absent, we can find that the P — V 
criticality of a five dimensional spherical Gauss-Bonnet black hole has an impressive analogy 
with the P—V criticality of a van der Waals liquid-gas system. This is a non-trivial result, 
since the previous studies without Gauss-Bonnet term in the literature indicate that there 
is no P — V criticality for the case of vanishing charge. Here the Gauss-Bonnet coefficient 
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Figure 2: The Gibbs free energy as a function of temperature for different pressures. The 
left plot is for the case with d = 5, where the “swallow tail” appears and a critical point 
exists. The right plot is for the case with d = 6, where there is no phase transition. In 
both plots we set a = 1. 


Q in some sense plays the role of charge. Clearly the existence of P — V criticality and 
the small black hole/large black hole phase transition is closely related to the existence of 
the small local stable Gauss-Bonnet black holes in AdS space [26]. Note that such a stable 
small black hole does not exist when d > 6. 

In this subsection, we will study the critical behavior of some physical quantities and 
compute corresponding critical exponents. As a warm up, we will first give a brief intro- 
duction to the critical exponents for the fluid/gas phase transition in the van der Waals 
liquid-gas system before applying it to the black hole system. 

Near the critical point, the critical behavior of a van der Waals liquid-gas system can 
be characterized by the following critical exponents [30]: 


C, ~ (=a 
wae „Toiy, 
Ve Te 
kr ~ E 
P-P. ~ wa (29) 


where T < T, and C, = T(23)|, is the the heat capacity at constant volume and Kp = 
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—v '(%) Ir is the isothermal compressibility. The index c denotes that a quantity is taken 
value at the critical point of the van der Waals liquid-gas system. Subscripts g and / stand 
for the gas phase and liquid phase, respectively. For the van der Waals liquid-gas system, 
the critical exponents take values as follows: a = 0, 8 = 1/2,y=1,6 = 3. 

Now let us turn to calculate the critical exponents in the Gauss-Bonnet black hole 
system. In order to examine the critical behavior near the critical point, we introduce the 
following expansion parameters 


E 


Then we can make Taylor series expansion for the equation of state (17) with d = 5 and 
k = 1. We assume the expansion has the form as [8] 


p = 1+ aioT + ai TE + ao + Ole, e$). (31) 


We can obtain all the expansion coefficients in (31) through computing its derivatives with 
respect to T and € at the critical point. These coefficients turn out to be 


ajo = 4, a4, = —6, a03 = —6. (32) 


Using Maxwell’s equal area law we obtain the following equation 


“s ap 3 1 
0= / cade => qos (€ — ef) + saute —¢)=0. (33) 
€l 


On the other hand, note that during the phase transition the pressure remains unchanged, 
we then have 
p'le = P'le = aT (€g — €1) + ao3(63 — €) = 0. (34) 


These two equations (33) and (34) will be used to calculate the critical exponents here 
and in the next section where the electric charge does not vanish. Now substituting these 
expansion coefficients in (32) into these two equations, we conclude that these equations 
have a unique non-trivial solution (€; 4 €,) only when 7 < 0. In other words, when T < Te, 
there are two different volumes for a same pressure. The nontrivial solution of (33) and 
(34) with (32) reads 

a&=-vV-T, €=V-T. (35) 


From (29) we can immediately conclude that the critical exponent 6 = 1/2. The isothermal 
compressibility can be calculated as follows. 


1 Ov 1 1 
Z — = 36 
i OP in. “= p 6T’ (Pe) 


which indicates that the critical exponent y = 1. In addition, it can be easily seen that 
(p’ —1)|,=0 = —6e?, which tells us 6 = 3. And the heat capacity near the critical point has 
C, = T(23),, = 0, we then have a = 0. 
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Thus we have obtained all the critical exponents in the P — V criticality of a five- 
dimensional spherically symmetric Gauss-Bonnet black hole in the extended phase space 
by treating the cosmological constant as a pressure. It is easy to check that these critical 
exponents satisfy the following thermodynamic scaling laws 


a+t26+y=2, a+ 6(14+6)=2, 
(1 +8) = (2-a) - 1), y= (6-1). (37) 
It has been confirmed that these critical exponents associated with the P — V criticality 


of the black hole are the same as those in the van der Waals liquid-gas system, and they 
obey the same scaling laws. 


4 P-V criticality in the case with Q # 0 


In this section, we will consider the case with non-vanishing charge Q, i.e., the charged 
Gauss-Bonnet black holes in AdS space. In this case, we have the equation of state from 
the temperature (7) of the black holes as 


d—2 2kha d—2)(d—3)k d — 2)(d —5)k?a Q 
p=- 7 ae eee. ey 
Th ri 16rri 16771, SIT), 


We see that when d = 5, the third term identically vanishes in the pressure. Therefore the 
behavior of the pressure in the case of d = 5 is different from the case with d > 5. Thus we 
will discuss these two cases separately. In addition, for later discussions, it is convenient 
to rescale some quantities in the following way 


Be se -2 =2 a1 
Th = 2/Q|*3, a = b|Q|7 2 P = p|Q|7 3; T = t|Q|7 3, G=4g/Q|, (39) 


where d is the dimension of the black holes and G the Gibbs free energy. 


4.1 Charged black holes in five dimensions 


When d = 5, the equation of state reduces to 


3 2ka 3k Q? 
P = — (1 + —> T - — + —. 40 
TER T if ) 8rr? T Snr? (40) 


The critical point is determined by the equation (20). The critical point has to satisfy the 
following equation 


rp, — 6kr,,.a — 5kr? Q? — 18aQ? = 0. (41) 
And the critical temperature can be written as 
k 4 _ N2 
NN (42) 


mrs (r2,+6ka) 
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Rewriting these equations in terms of dimensionless quantities in (39), they become 


3 j 2h 3k 1 


= — a eee 4 
P ie x? ) 8ra?  8rrê’ (43) 
xÊ — 6kba, — Ska? — 18b = 0, (44) 
kx? = 1 
ae (45) 


max3(x2 + 6kd) | 
It is easy to see that for the Ricci flat black hole with k = 0, the pressure is a monotonic 
function of the horizon radius x. There does not exist any phase transition and therefore no 
critical point exists in this case. In the following we will focus on the cases with spherical 
topology (k = 1) and hyperbolic topology k = —1, respectively. 


1) spherical case k=1 
In this case, the equation (44) reduces to 
xê — 6brt — 5x? — 18b = 0. (46) 


Before solving this equation exactly, let us take a little time to consider the polynomial in 
the left hand side of (46). We can easily get some qualitatively useful information. When 
b > 0 (which means @ > 0), the equation admits only one positive real root denoted by 
£e. Furthermore, since the value of the polynomial at £e = 1(which means r, = |Q]) is 
—4 — 24b which is definitely negative for b > 0, so we can conclude that the positive real 
root of (46) must satisfy xe > 1, which leads to a positive critical temperature (see (45)). 
Although this equation can have an analytic expression for the solution, we will not present 
it here. The p — x diagram and Gibbs free energy are shown in figure (3). 

In the p — x diagram, the critical isotherm T = T, is denoted by the red line, where a 
second order phase transition happens just as the case of the van der Waals system. The 
curves above the critical line have a positive compression coefficient, corresponding to the 
“gas” phase. Each curve below the critical one corresponds to the case having temperature 
lower than T}. Such curves have two stable branches, one is for small x and the other is for 
large x, and the middle one between them is unstable. At these temperatures, the small 
black hole/large black hole phase transition occurs, which can be seen more clearly in the 
right plot. The Gibbs free energy as a function of temperature develops a “swallow tail” 
for P < P., indicating a typical first order phase transition. The situation discussed here 
is comparable to the case without charge. We therefore expect that this phase transition 
could have the same critical exponents as the case in Q = 0, which is indeed true as shown 
shortly. 


2) hyperbolic case k=-1 
When k = —1, the equation (44) reduces to 
xÊ + 6ba, + 5x? — 18b = 0. (47) 
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Figure 3: The p— «x diagram (left plot) and Gibbs free energy (right plot) near the critical 
temperature and critical pressure for the case of k = 1 in five dimensions. In this figure we 
take b = 0.2. pe and te are obtained by Eqs. (43),(44), and (45) with k = 1 and b = 0.2. 


This equation has at most one real positive root denoted by xe. It can be seen that this 
root must be located in the range (0,3,/2b/5). Note that when k = —1, the non-negative 
definiteness of black hole entropy (8) gives the constraint on the horizon radius, 


x, 2a=4 


rj, + =P ha > 0. (48) 


As d= 5 and k = —1, the above constraint leads to 
x? > 6b. (49) 


With this constraint, one can see that the critical temperature (45) is always negative. 
Thus we conclude that in the case k = —1, there does not exist any critical point and no 
phase transition happens. We can also obtain this conclusion from (43) by noting the fact 
that when k = —1, the pressure is a monotonic function of horizon radius. 


4.2 Charged black holes in more than five dimensions 


In this case, d > 6, the equation of state of the black holes is give by (38). The critical 
point has to satisfy the equation 


(d —3)a2*-* — 12kba?4-* + 12(d — 5)b?a2*-® — (4d — 10)kz? —12(2d—7)b=0. (50) 
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And the critical temperature is given by 


pa k(d — 3)x?®™ ° + (2d — 10)ba24-8 — 2 i 

2nx24-7 (x2 + 6kb) on 
We can see from (38) that when k = 0, the pressure is always a monotonic function of the 
horizon radius for a fixed temperature. Therefore in the case of k = 0, the “gas” (large 
black hole) phase always dominates and it is impossible to have critical point and phase 
transition. We will discuss below the cases of k = +1 separately. 


1) spherical case k=1 


When k = 1, one can see that the equation (50) has only one positive root. The critical 
point, if any, has to satisfy the following constraints, 


(d — 3)a2¢-* — 12ba74-* + 12(d — 5)b?24-® — (4d — 10)x? — 12(2d — 7)b = 0, (52) 
(d — 3)a24-6 — 2 + (2d — 10)bz?4* > 0, (53) 
F|T=T r= = 0. (55) 


The constraint (53) comes from the requirement to have a positive critical temperature. 
And the constraint (54) comes from (5), which insures the black hole solution to be asymp- 
totical AdS. In the case of Q = 0, we see that the Gauss-Bonnet term plays the role of the 
charge in P—r, diagram. When b is very large (b = a|Q|~2/“4-9) >> 1), which means @ will 
be much large than Q, one may expect that in this case, the situation is the same as the 
case of Q = 0, where there does not exist a physical critical point when d > 5. This can 
be seen from the pressure (38). So we conclude that there could exist a physical critical 
point only when b is not too large. This argument can also be supported by analyzing 
(52). When b is very large, one has an approximate solution of the equation (52) as 


~ (a, (56) 


while the pressure (38) is dominated by the third term as 


(d — 2)(d —5)b 


= 16rzź 


(57) 
Clearly in this case the pressure is negative and the constraint (54) cannot be satisfied. The 
exact range of b, in which a physical critical point exists, can be determined numerically. 
We plot the range of the parameter b in the cases with different dimensions in figure (4). 
We see from figure (4) that when the dimension of spacetime increases, the upper bound 
of b decreases, but finally it arrives at a finite one around 0.27. As a result, in the special 
case with b = 0, a critical point always exists when d > 6. This is in agreement with 
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k=1,|Q|>0 


Figure 4: The range of the parameter b in the cases with different dimensions. In this 
range, a physical critical point exists. 


the one for the charged AdS black holes analyzed in Ref.[9]. Clearly the existence of the 
upper bound for the parameter b is mainly due to the competition between the third term 
and the fourth term of the pressure (38) in the region of small ra. When the third term 
dominates, the critical point is absent, while it will appear if the fourth term dominates. 
Therefore the existence of the upper bound of b is the effect of the Gauss-Bonnet term. 


2) hyperbolic case k=-1 
When k = —1, the physical critical point should satisfy the following constraints, 


(d — 3)a24-4 + 12ba?4- 4 12(d — 5) 6224-8 + (4d — 10)x? —12(2d-—7)b=0, (58) 


[(d — 3)x24-6 + 2 — (2d — 10)br?™’] (x? — 6b) < 0, (59) 
F| < 0, (61) 
>». 2d—4 4 
> = —— )b. 
t 2 (2+ 777) (62) 


The former four constraints have the same origins as in the case of k = 1. The fifth one 
(62) comes from the requirement that the black hole entropy must be non-negative. It is 
easy to see that the equation (58) has only one positive root. Before discussing the generic 
case, let us discuss two extremal cases. 
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When b is very small, we can find that the solution of (58) is x. ~ \f SEP. Both 


constraints (59) and (62) can be satisfied. But the constraint (60) cannot be satisfied, 


because in this case, 
b 1 


~ Sa ~ 8rbt-3 
as b + 0. Therefore we can conclude that when b is small enough, there are no phase 
transition and critical point. This is consistent with the observation that when b = 0, 
there does not exist any real root of equation (58). 

On the other hand, when b is very large, one can have an approximate solution of (58) 


as r 
Qd—7 == 


bpe => œ, (63) 


b(d — 5) 


In that case, we can see that the constraint (62) cannot be satisfied. 

As a result, one may expect that those constraints can be satisfied only when 6 is in 
some region. However, we can numerically check that the real positive root of equation (58) 
is out of the region to have a positive Hawking temperature of the black holes. Therefore 
we conclude that the critical point is also absent in this case. This is consistent with the 
fact that the heat capacity of the black holes at a fixed charge is always positive in this 


Case: 
— 2), rt (r2 + 2ka)T (OT \ 
Cee (Sr) _ (d —2)der,, wi + 2kã) (=) (64) 
P,Q Th? PQ 


where (£) is always positive in the physical region of the horizon radius. 
P,Q 


4.3 The critical behavior when d > 5 


We can calculate the critical exponents of some physical quantities for Q Æ 0 in d > 5 in 
the same way as the case for the black holes without charge in d = 5. It is found that 
those critical exponents are universal and does not change. That is to say, those critical 
exponents in the case of charged black holes in d > 5 are the exactly same as those in the 
case of natural black holes in d = 5 calculated in Section 3. We therefore will not repeat 
these calculations here, instead we pay some attention on some special aspects. Near a 
critical point, the pressure behaves as (31). From the requirement that the pressure should 
be equal at transition point, we have 


diiT (€g — €1) + TAG -&)=0. (65) 


On the other hand, from the Maxwell’s equal area law, one has 


3 
ay T(€ =E) F 5203 (€5 =g) =0. (66) 
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Then nontrivial solution of the above two equations appears only when a1,a93T < 0, which 


is given by 
y —411403T y —411403T 


g= €=— 
3] aos 3] aos 


(67) 
In section 3, we have obtained the exact values of a1; and ag3 in the case Q = 0 and d = 5. 
In a general case, it seems impossible to give an analytical expression for aj, and aog. As 
a result, we will numerically give two examples to show their behaviors. Figure (5) shows 
the case of d= 5, k = 1 and Q 4 0. We can see from the figure that both aj, and aoz are 
negative. Thus we can conclude that the phase transition can occur when the temperature 
is less than T}. In addition, let us notice that when b increases, both a, and ao3 quickly 
approach to —6, the value in the case of Q = 0. This is consistent with the fact that when 
b increases, the effect of Q decays quickly by noting the relation |Q| = (a/b) in d = 5 
dimensions. 


d=5,k=1 |Q>-0 


Figure 5: The values of a1, and aog when d = 5, k = 1 and Q 4 0. They are both negative. 
So the phase transition can occur when T < 0. 


Figure (6) shows the case of d = 7, k = 1 and Q # 0. We can see that in this case, 
the behaviors of a1; and dog are different from the case of d = 5. Both a4, and aoz do 
not approach to a fixed value. When b is less than its upper bound, the critical point is a 
physical one, in this case the phase transition can occur when T < Te since ay, < 0 and 
ao3 < 0; when b is beyond the upper bound, there do not exist any critical point and phase 
transition. 


5 Conclusions and discussions 


In this paper we studied the phase transition and critical behavior of topological charged 
Gauss-Bonnet black holes in d-dimensional AdS space. We discussed this issue in the 
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d=7,k=1,|Q|>0 


Figure 6: The values of aıı and aoz in the case of d = 7, k = 1 and Q # 0. In the region 
where a1, and aoz are both negative, the phase transition can occur when T < 0. 


extended phase space where the cosmological constant is treated as the pressure of the 
thermodynamic system and its conjugate quantity is the thermodynamic volume of the 
black holes. The black hole horizon can be Ricci flat (k = 0), spherical (k = 1), or 
hyperbolic (k = —1), respectively. By scaling argument, we obtained a generalized Smarr 
relation (13) for the Gauss-Bonnet black holes. This relation implies that one has to 
introduce the conjugate quantity A to the Gauss-Bonnet coefficient & in the first law of 
black hole thermodynamics (11). We calculated the conjugate quantity A and found that 
it consists of two terms, one term comes from the contribution of the Gauss-Bonnet term 
to the mass of the black holes, while the other from the contribution of the Gauss-Bonnet 
term to the entropy of the black holes. That is to say, one term is related to the volume 
of the black holes and the other is related to the area of the black holes. 

We showed that for the Gauss-Bonnet black holes with a Ricci flat or hyperbolic horizon, 
no phase transition can appear and the large black hole phase (gas-like phase) is always 
dominated. This conclusion is related to the fact that for the Gauss-Bonnet black holes 
with Ricci flat or hyperbolic horizon, their heat capacity is always positive and therefore 
those black holes are always local thermodynamically stable [26]. In particular, for the 
Ricci flat case, those thermodynamic quantities of the charged Gauss-Bonnet black holes 
are completely the same as those of the black holes without the Gauss-Bonnet term [see 
(6)(7)(8)]. As a result, there also do not exist P — V criticality and small black hole/large 
black hole phase transition in the charged AdS black holes with Ricci flat or hyperbolic 
horizon in general relativity. Furthermore we believe that this conclusion is also valid for 
AdS black holes with Ricci flat or hyperbolic horizon in other gravity theories. 

For the Gauss-Bonnet black holes with a spherical horizon, we found that even when 
the electric charge is absent, the P—V criticality and the small black hole/large black hole 
phase transition can appear. It happens only in the case of d = 5, but not for the case of 
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d > 6. This is related to the fact that in the d = 5 case, there exists a local stable small 
black hole, while it is absent in the case of d > 6 [26]. Also we can see from the equation of 
state (17) that when d > 6, the last term always dominates in the region of small horizon 
radius, which implies that there is no branch with stable small black holes and that no 
phase transition can occur. When the electric charge is not vanishing, the P — V criticality 
and the small black hole/large black hole phase transition always appear in the case of 
d= 5. Namely in this case, there is no any limitation on the parameter b = 4|Q|?/(4-%). 
While d > 6, there is an upper bound of the parameter b, beyond which there do not exist 
any P — V criticality and phase transition. The upper bound is dimensionally dependent 
and we obtained numerically the upper bound in some dimension cases. This difference 
between d = 5 and d > 6 can be clearly understood from the equation of state (38). In 
the case of d = 5, the third term identically vanishes and thus last term related to the 
charge of the black holes always dominates in the region of small horizon radius. Note that 
the last term is positive and this implies that there exists a phase with stable small black 
holes. Therefore in this case the small black hole/large black hole phase transition always 
can happen. While in the case of d > 6, the third term no long vanishes and thus there is 
a competition between the third term and last term in the region of small horizon radius. 
If the third term wins, there is no phase transition because the third term is negative and 
there does not exist the branch with stable small black holes in this case. On the other 
hand, if the last term wins, the situation behaves very like the case of d = 5. Namely in 
this case, there exists the branch with stable small black holes. Thus the P — V criticality 
and the small black hole/large black hole phase transition can appear. This is just the 
reason why when d > 6, there is the upper bound on the parameter b. Clearly the bound 
on the parameter b is the effect of the Gauss-Bonnet term. 

In all the cases when the P—V criticality appears, we calculated some critical exponents 
and found that those critical exponents are the same as those in the van der Waals liquid- 
gas system. This further provides the supports for the analogy between the black holes in 
AdS space and the van der Waals liquid-gas system. 

Finally let us mention that when the Gauss-Bonnet term is present, the effective AdS 
radius is not the one / in (1), but [26] 


l 4a 

lef = V2 1+ 1— P’ 
This can be seen by considering the vacuum solution of the theory, namely the solution 
(4) with M = Q = 0. In this case, one may wonder whether one should define an effective 
pressure according to the effective AdS radius (68) through the relation (1) in discussing the 
thermodynamics of the Gauss-Bonnet black holes in the extended phase space. It turns out 
not true. This reason is that to derive the first law (11) of black hole thermodynamics and 
the generalized Smarr relation (13) in the Hamiltonian formulism where the thermodynamic 
volume is included as a conjugate quantity of the cosmological constant, the cosmological 
constant A comes into the process, rather than the effective cosmological constant. See, 
for example, [19]. 


(68) 
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